Two-parameter scaling of correlation functions near continuous phase transitions 
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We discuss the order parameter correlation function in the vicinity of continuous phase transitions 
using a two-parameter scaling form G(k) = k^ 2 g(k^,k/k c ), where k is the wave-vector, £ is the 
correlation length, and the interaction-dependent non-universal momentum scale k c remains finite 
at the critical fixed point. The correlation function describes the entire critical regime and captures 
the classical to critical crossover. One-parameter scaling is recovered only in the limit k/k c — » 0. 
We present an approximate calculation of g(x,y) for the Ising universality class using the functional 
renormalization group. 
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Correlation functions in the vicinity of continuous 
phase transitions usually assume a scaling form at long 
wavelengths*^ Asymptotically, for wave-vectors k much 
smaller than the relevant microscopic scale A (such 
as the inverse lattice spacing), the order parameter 
correlation function G(k) can be written as G(k) — 
k~ 2+v g ± (fc£) , where the anomalous dimension 77 charac- 
terizes the universality class of the system, £ is the order 
parameter correlation length, and the scaling functions 
g + {x) and g~(x) describe the regime above and below 
the critical temperature T c , respectively. Precisely at 
the temperature T — T c there should be a unique scaling 
function, so that g + (oo) = g~(oo). If the dimensionality 
D of the system is smaller than its upper critical dimen- 
sion, then g (x) approach finite limits for large Ao£ and 
we may take the limit Ao — > 00. According to the single- 
parameter scaling hypothesis, the functions g^{x) are es- 
sentially determined by their asymptotic limits for small 
and large x, so that an extrapolation to the crossover 
regime x « 1 is possible either from the region x 3> 1, or 
from a: C 1. However, it is clear that another prominent 
scale, which is deducible already from a dimensional anal- 
ysis of Ginzburg-Landau-type theories but is absent from 
the one-parameter-scaling picture, will be important at 
larger wavevectors: the interaction-dependent scale k c 
which measures the size of the Ginzburg critical region. 3 
Here we shall assume a microscopic model such that the 
scale k c is small compared to the natural cutoff of the 
model. This happens e. g. for Ising models if the range of 
interactions is made very large^ or near the critical point 
of complex fluids. 7 In these systems a classical-to-critical 
crossover can be observed from a region dominated by 
the Gaussian fixed point to the region governed by the 
Wilson-Fisher fixed point. While the crossover of static 
thermodynamic derivates have been well studied in the 
past^^i 7 . the crossover behavior of the order parameter 
correlation function has only been investigated within a 
one-parameter theory valid right at T C &&1Q 

In this work we extend the one-parameter scaling the- 
ory for correlation functions to account also for the 
interaction-dependent scale k c <C Ao which divides the 
critical regime <C k <C Aq into two separate regimes 
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FIG. 1: (Color online) Macroscopic domain of wave- vector 
k -C Ao and correlation length £ _1 <C Ao. In the four shaded 
regions the correlation function G(k) has different behavior. 
The critical regime £~ <C k <C Ao is subdivided into a critical 
long-wavelength (CL) regime £ <C k <C k c (where the usual 
one-parameter scaling with anomalous dimension r\ is valid) 
and a critical short-wavelength (CS) regime <C k c <C k 
(where the self-energy exhibits universal power-law scaling 
distinct from the CL regime). This figure is a refinement of 
the corresponding Fig. 1 of Ref. 



where G{k) has rather different properties ; 8 ' 9 ' 10 only in 
the critical long-wavelength (CL) regime <C k -C 
k c <C Ao does the correlation function scale asymp- 
totically with an anomalous dimension rj. For /c c £ 3> 
1 there exists another critical short- wavelength (CS) 
regime C fc c < C Ao where the behavior of G(k) 
is still universal but distinct from the anomalous scaling 
in the CL regime. The scale k c is present in the usual 
renormalization group (RG) analysis^ but is lost if the 
RG flow is linearized in the vicinity of the critical fixed 
point. 

Taking into account the macroscopically ordered and 
disordered regimes, four different macroscopic domains 
of wave-vectors should be distinguished, as summarized 
in Fig. [I] In order to bring out the difference between the 
CL and the CS regime, we express G(k) in terms of the 
irreducible self-energy £(fe) = G^ 1 (k) — k 2 . In the CL 



regime scales asymptotically as k 2 ~ v which domi- 

nates the bare fc 2 -dispersion, so that G(k) ~ 
On the other hand, in the CS regime S(fc) may or may 
not be larger than fc 2 , depending on the strength of the 
interaction. What is more important here is that for 
k c <C k -C Ao the self-energy can still be expressed in 
terms of a universal scaling function which shows power 
law behavior different from the CL regime. Therefore, 
the self-energy in the macroscopic domain k <^ Aq should 
be written in terms of scaling functions cr ± (a;, y) depend- 
ing on two parameters x — fc£ and y = k/k c , 



E(jfc) = k^ik^k/K) 



(1) 



where o~ + (x,y) corresponds to the disordered phase T > 
T c , while a~(x,y) describes the ordered phase T < T c . 
If the system is cooled through T c and one probes the 
system at a fixed scale k, then the asymptotic critical be- 
havior with anomalous exponent r\ can only be observed 
in the CL regime <C k <C k c . Yet, the same exper- 
iment at k c <C k would reveal the universal behavior of 
the self-energy in the CS regime. The non-universal scale 
k c , defining the width of the CL regime, is missed within 
the field-theoretical renormalization group, which effec- 
tively sets k c — oo in the self-energy. However, k c and 
the scaling functions a ± (x, y) can be calculated using the 
functional renormalization groupe d 2 ' 13 ! 14 ' 15 

We believe that the statements above are general and 
apply to any continuous phase transition. In the rest 
of this work we shall demonstrate their validity by an 
approximate calculation of the scaling functions a (x, y) 
for the Ising universality class in D dimensions, which 
can be modeled by an action involving a real field <p(r), 
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where an ultraviolet (UV) cutoff Ao is assumed to regu- 
larize the theory. It is instructive to consider first the per- 
turbative calculation of the self-energy, which is possible 
as long as the relevant dimensionlcss coupling constant 
Uo = ua £ 4 ~-° is small. The correlation length £ can be 
expressed in terms of the self-energy as £~ 2 = £>E(0), 
where the field renormalization factor Z is defined as 
Z- 1 = 1 + dE(k)/dk 2 \ k=0 - For D < 4, i.e. below 
the upper critical dimension of our model, the loop in- 
tegrals generated in the perturbative expansion are UV 
convergent and we may take the limit Ao — > oo. Per- 
turbation theory yields the self-energy in scaling form, 
£ 2 E(/c) = Z- 1 + Aa^(k£). The lowest order dia grams 
giving rise to a momentum dependence of the self-energy 
are shown in Fig. [21 In the disordered phase the corre- 
sponding scaling function is (up to order 0(uq)) 



A«t+(x) = ^ / X (P) 
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where n is an arbitrary unit vector, J p = J ; and 
X(p) = j p , {[p' 2 + iW+pf + I]) -1 - For T<T C there 




ln(x) = In (k£) 

FIG. 2: (Color online) Lowest order perturbative results for 
the scaling functions Aa^ (x) (solid line) and Ao~g (x) (dashed 
line) in D = 3, see Eqs. © and ©. The function AcrJ (Actq ) 
is determined by the upper (lower) diagram on the left. 
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is a finite three-legged vertex of order u , so that the 
leading momentum dependent contribution to the scaling 
function A(7q~ (x) is given by the lower diagram in Fig. [2] 
which is linear in Uq, 



Aa (x) 



3u 



[ X (0)-X(x)]+O(u 2 ). 



(4) 



An explicit evaluation of Acr^ (x) in D = 3 is shown in 
Fig. O The qualitative behavior of the functions cr ± (x) 
is easily obtained for arbitrary D. While Aa^(x) oc x 2 
for x — > 0, the asymptote for large x is non-trivial due to 
the non-analytic behavior of the function x(l) ^ q D A 
for large q in D < 4. We find Aa^(x) cx X 2 ^ D ~ 3 ^ for 
3 < D < 4, Acr f (x) oc hxx for D = 3, and Aa^(x) ~ 
A<Jq (oo) = 0(1) for 2 < D < 3. In the ordered phase, 
the scaling function Aa^ (x) approaches for large x a fi- 
nite limit Ao-q (oo) = |uox(0) but the sub-leading cor- 
rection is non-analytic, Aa^ (x) — A(7q~(oo) oc x d ~ 4 . 

Lowest order perturbation theory does not directly 
reveal the existence of a second characteristic length 
scale fe" 1 besides £. However, the perturbative approach 
breaks down in the vicinity of the critical point for D < 4 
since the effective dimensionless expansion parameter 
i*o = ua ii ( 4_i) diverges. A consistent resummation of 
these divergences using a RG approach not only leads 
to non-trivial renormalization of £, but also reveals the 
presence of the characteristic scale k c £ The signature of 
both scales k c and £ -1 is already visible in the usual 
one-loop RG equations for the rescaled coupling param- 
eters Ti = Zita/A 2 and ui = KjjZfA D ^ i UA, describing 
the evolution of r\ and ma in Eq. (J5]) as the degrees of 
freedom in the momentum shell A = Aoe~ ( < k < Ao 
are integrated outt 1 "" 2 - Here, Zi is a field renormalization 
factor and K D = 2 1 - d tt- d / 2 /T(D/2). For T > T c the 
one-loop RG equations are 1 - 

dm = 2n + m/[2(i + rj)] , (5) 
d m = (4-£) U/ -3M 2 /[2(i + r ; ) 2 ]. (6) 

If we fine-tune the initial value ro such that the system is 
close to the critical point, the typical RG flow of r; and 
ui as a function of the "RG time" I is shown in Fig. [3] 
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FIG. 3: (Color online) Flow of the coupling parameters ri and 
ui for a nearly critical system as a function of the RG time I 
obtained from Eqs. ([5]) and (|6} for D — 3. 



The two characteristic RG times l c and Z* clearly separate 
three distinct regimes. In the regime < Z < Z c , the cou- 
pling parameters r; and ui flow toward the values r* and 

associated with the critical fixed point; in the interme- 
diate interval l c < I < Z* the RG flow is close to the fixed 
point and very slow while for Z* < I the trajectory rapidly 
flows away from the fixed point. A simple calculation 9 
yields the estimate l c « (4 — D)^ 1 ln(u*/uo) for uo <C u*. 
The one-parameter scaling hypothesis is based on the 
analysis of the linearised RG flow in the vicinity of the 
fixed point and correlates with the existence of only one 
unstable direction. The eigenvalue corresponding to the 
unstable direction then directly determines the scaling 
of £, see e. g. Refs. [HQ- The underlying assumption is 
that as long as the initial parameters are close to the 
critical surface (not necessarily close to the fixed point), 
the unstable direction will completely dominate the flow 
and thus the behavior at small k. Corrections to pure 
scaling are known to arise, if the stable (i.e. irrelavant) 
directions of the linearised flow are also accounted for^iiii 
However, the analysis leading to one-parameter scaling is 
incomplete since the flow on the critical surface itself will 
typically introduce (at least) one additional scale which 
is not visible in the linearized flow. This is how the char- 
acteristic size k c of the Ginzburg critical region enters 
and how the one-parameter scaling hypothesis can be 
extended*- Below, we show that the RG times l c and Z* 
and the associated scales k c = Aoe~' c and = Aoe~ l * 
separate regimes where the momentum dependent self- 
energy shows qualitatively different behavior. 

To obtain the self-energy in the vicinity of the crit- 
ical point taking into account the full RG flow we use 
the functional renormalization group (FRG)^ ' 12 ' 13 ' 14 ' 15 
An exact hierarchy of flow equations for the one-line ir- 
reducible vertices of our model is obtained by differen- 
tiating the corresponding generating functional with re- 
spect to A and then expanding this functional in powers 
of the fields. 13 ! 15 The true self-energy is then obtained 
as E(fc) = limA^o ^a(/c). To calculate the scaling func- 
tions defined in Eq. ([I]), it is convenient to work with 
rescaled variables which manifestly exhibit the scaling 
dimensions of all parameters. With dimensionless mo- 
menta q = k/A we define the dimensionless self-energy 




FIG. 4: (Color online) Graphical representation of Ti(q) in 
Eq. |[7J describing the effect of higher order irreducible ver- 
tices, represented by numbered circles, on the evolution of 
the self-energy as degrees of freedom are integrated out. The 
crossed circle represents the flowing order parameter, small 
black dots represent the derivative with respect to A, circles 
with G [G] the exact [single scale] propagator, and Si, 2 sym- 
metrizes with respect to external labels. For T > T c , the last 
two diagrams involving the three-legged vertex vanish. 



as Ti(q) = ZiA^ 2 T, A (Aq), where I = -ln(A/A ). It sat- 
isfies an exact flow equation of the form 9 - 



d l r l (q) = (2- m -qd q )T l (q)+t l (q) , 



(7) 



where r\i 



-dilnZi — dti(q)/dq 2 \ q= o is the flowing 



anomalous dimension. The function Ti(q) depends on 
higher order irreducible vertices and describes the usual 
mode elimination step of the RG procedure. A graphi- 
cal representation of the interaction processes contained 
in Ti(q) is shown in Fig. 2J Note that in the ordered 
phase, where our field has a finite vacuum expectation 
value (<p k ) = (2n) D S(k)M A , Eq. has to be aug- 
mented by a flow equation for the flowing order param- 
eter Ma, which is obtained by requiring that the ver- 
tex with one external leg vanishes identically for all AJ£ 
Defining 7/(g) = Ti(q) — T/(0), the physical self-energy 
can be written as an integral over the entire RG trajec- 
tory, 



= hm 
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dUe-v+ft^Tt&k/Ao). (8) 



To relate this to the scaling functions cr ± (a;, y) defined in 
Eq. ([1]), we use the definitions k c = Aoe _ic , x = y 
k/k c . By construction e~ 21 * = (Ao£) -2 = lim; 



— 2/ 

e u n 



and x/y = (fc c £) 



J, -la 



Our final result is 



cr ± (x,?/) 



JT 

Zx 2 



+ / dle-W-U+fi^itf-^y). (9) 



The corresponding scaling form of the order parame- 
ter correlation function is G(k) — k~ 2 g ± (k£,k/k c ), with 
g ± (x,y) = [y 2 + cr ± (a;, y)] _1 . The behavior of g ± (x 1 y) 
crucially depends of the ratio y/x. Exactly at the crit- 
ical point one has x = 00, so that <J*(y) = «t ± (cxd, y). 
In the CL regime y -C 1 we obtain the usual anoma- 
lous scaling <r*(?/) oc y 2 ~ r ', while in the CS regime y ^S> 1 
<J*{y) oc lny in D — 3, which resembles the perturbative 
result (u* / uq) 2 A<Jq (k / k c ) shown in Fig. [2] where k c now 
plays the role of an effective infrared cutoff^ 

Away from criticality, Eq. §§§ can only be evaluated 
numerically. The result of a calculation of a + (x,y) for 
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FIG. 5: Scaling function Aa + (x,y) = a+(x,y) - (fc^)" 2 ^" 1 
on the line x/y — k c £ w 1700 in D = 3 using a truncation 
of the FRG flow equations similar to the one used in Ref . [|| . 
Also shown are (tt*/uo) 2 Aoo + '(fc/fc c ), which gives a reason- 
able approximation of the CS regime, and the asymptotic crit- 
ical behavior y 2 ~ v (our truncation leads to r/ ~ 0.101) which 
is approached within the CL regime before it is replaced by a 
y 2 behavior in the macroscopic disordered regime. 

our model in D — 3, using a similar truncation of the 
FRG flow equations as in Ref. 9 , is shown in Fig. [5l the 
crossover between the different regimes is clearly visible. 
A FRG calculation of <r~(x, y), using the flow equations 
for the irreducible vertices in the ordered phase derived in 
Ref. [15j |. will be given elsewhere^ In the regime j/>i, 
which corresponds to k c <C the system is far away 
from criticality (see Fig. [T]), and the scale k c disappears 
from the problem. In this case it is better to normalize 
the self-energy differently, 

/>oo 

fE(fc) = Z- 1 + / dle- 2 V- l ^+f° dTr i-ji(e l - l *x). (10) 
Jo 



Evaluating this expression to lowest order in the inter- 
action using the trivial scaling r; rs e 2 'r; = o and ui ~ 
e (4-- D )' U j =0! the second term on the right-hand side re- 
duces for T > T c to the function Acr^ (x) given in Eq. © 
and for T < T c to the function A<Jq (x) given in Eq. ([4|. 

In summary, we argue that the correlation function 
G(k) near continuous phase transitions has a scaling 
form which involves (at least) two parameters. While 
in the critical long wavelength (CL) regime the self- 
energy shows asymptotically the usual anomalous scal- 
ing, S(fc) oc A; 2-1 ', there exists another macroscopic crit- 
ical short-wavelength (CS) regime k c -C k <C Ao where 
the scaling of the self-energy is given by a different power 
law. The scaling function describes the crossover be- 
tween these regimes and the macroscopically ordered and 
disordered regime, as well as the finite-fc corrections to 
the asymptotic k 2 ~ v behavior (which is reached only for 
k — > 0) in the CL regime. The non-universal scale k c 
derives from the nonlinear behavior of the RG flow on 
the critical surface and is missed if the flow is linearized 
around the critical fixed point. Although for concreteness 
we have focused here on the Ising universality class, the 
two parameter scaling applies to the critical behavior of 
systems which in the continuum limit are well described 
by a Ginzburg-Landau-type action in which terms which 
are higher order than quartic in the fields are negligible. 
The CS regime k c <C k <C Ao is accessible only if k c is 
small compared with the effective UV cutoff Ao. For the 
case of generalized Ising-models, this requires that the 
range of interacions is very large A similar crossover 
might be more readily observed in near-critical polymer 
solutions where the crossover scale depends on the radius 
of gyration of the polymer^ which can be very large. 

We thank E. Vicari and A. Pelissetto for interesting 
and useful correspondence. 
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